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Abstract 

In  this  paper,  the  representation  of  the  kinematics  of  inelastic  flow  problems  involv¬ 
ing  finite  strains  and  rotations  is  discussed.  The  concept  of  multiplicatively  decompos¬ 
ing  the  deformation  gradient  into  elastic  and  plastic  components  is  utilized,  where  the 
plastic  deformation  gradient  represents  a  stress-free,  relaxed  configuration.  Here,  it  is 
demonstrated  that  the  choice  of  relaxed  configuration  is  not  essential  in  the  problem 
solution.  Therefore,  the  kinematic  decomposition  of  elastic-plastic  deformation  may  be 
chosen  to  best  analyze  the  specific  material  model  of  concern.  This  is  demonstrated  by 
analyzing  two  specific  materials,  the  planar  single  crystal  and  the  glassy  polymer,  using 
various  kinematic  representations. 

1  Introduction 

The  representation  of  inelastic  material  behavior  under  finite  strains  and  rotations  has 
been  the  subject  of  much  discussion  in  the  recent  plasticity  literature  [FARDSHISHEH 
AND  ONAT  (1972),  ONAT  (1982),  LORET  (1983),  DAFALIAS  (1984,  1985),  ANAND 
(1985),  AGAH-TEHRANI,  ET  AL.  (1987)].  Finite  strain  formulations  generally  be¬ 
gin  with  a  description  of  the  deformed  body.  This  may  be  given  by  the  deformation 
gradient,  F  =  Vxx  where  X  represents  the  reference  position  of  a  material  point,  x 
represents  the  current  position,  and  V  is  the  gradient.  The  multiplicative  decomposi¬ 
tion  of  the  deformation  gradient  into  elastic  and  plastic  components  was  first  proposed 
by  LEE  (1969),  thereby  introducing  the  concept  of  a  relaxed,  intermediate  configuration 
represented  by  the  plastic  deformation  gradient.  Recently,  the  physical  significance  of 
such  an  intermediate  configuration  has  been  debated.  HILL  AND  RICE  (1972),  ASARO 
(1983),  HAVNER  (1982)  have  proposed  a  convenient  framework  for  the  representation 
of  single  crystal  plasticity  where  the  rotation  of  the  crystal  lattice  is  perceived  to  be  an 
elastic  deformation.  The  plastic  deformation  gradient  represents  the  relaxed  configu¬ 
ration  such  that  the  crystal  is  stress-free  and  the  lattice  is  in  its  reference  orientation. 
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LORBT  (1983)  and  ANAND  (1985)  have  extended  this  basic  framework  for  the  single 
crystal  to  polycrystalline  materials  by  the  utilization  of  the  concept  of  Mandel  director 
vectors  to  monitor  material  orientation  and  the  application  of  representation  theorems. 
Here,  the  relaxed  configuration  represents  the  configuration  obtained  by  elastically  un¬ 
loading  such  that  the  director  vectors  are  in  their  reference  orientation.  In  this  case, 
one  faces  the  difficulty  of  assigning  a  set  of  director  vectors  to  each  material  element, 
whereas,  for  the  single  crystal,  the  vectors  are  naturally  defined  by  the  lattice.  LEE 
(1969)  and  FARDSHISHEH  AND  ONAT  (1972)  pointed  out  that  a  unique  relaxed  config¬ 
uration  can  be  obtained  by  elastically  unloading  to  a  stress  free  state  without  rotation, 
in  the  polar  decomposition  sense.  This  second  form  of  decomposition  has  the  advan¬ 
tage  of  being  quite  general,  but  it  does  not  explicitly  include  the  orientation  effects 
that  occur  in  crystalline  materials  and  which  are  included  in  the  framework  of  the  first 
decomposition,  where  the  relaxed  configuration  is  determined  by  the  material  orienta¬ 
tion.  However,  ONAT  (1982)  has  theoretically  shown  that  the  second  representation 
can  accomodate  the  internal  state  and  orientation  of  a  material  element  by  the  up¬ 
dating  of  appropriate  scalar  and  tensor  state  variables,  thereby  bypassing  the  need  for 
director  vectors. 

In  this  paper,  it  is  demonstrated  that  the  problem  of  the  single  crystal,  which  has 
been  examined  utilizing  the  first  decomposition  [PIERCE,  ASARO,  AND  NEEDLEMAN 
(1983)],  can  also  be  examined  utilizing  the  second  decomposition  together  with  the 
monitoring  of  the  orientation  of  state  with  tensor  state  variables  as  suggested  by  Onat. 
The  representation  of  the  inelastic  deformation  of  glassy  polymers  is  also  examined.  In 
this  analysis,  two  separate  but  unique  kinematical  decompositions  and  the  constitutive 
model  proposed  in  PARKS,  ARGON,  AND  BAGEPALLI  (1984)  and  BOYCE,  PARKS, 
AND  ARGON  (1988)  are  employed.  These  examples  of  the  single  crystal  and  the  glassy 
polymer  utilising  a  variety  of  kinematic  representations  demonstrate  that  a  particular 
choice  of  the  relaxed  configuration  is  not  essential  for  the  problem  solution.  However, 
choosing  a  convenient  relaxed  configuration,  specific  to  the  material  of  concern,  can 
simplify  the  solution  process  as  well  as  provide  a  physical  interpretation  of  the  deform¬ 
ing  solid.  As  will  be  shown  below,  the  "best”  choice  of  kinematic  representation  will 
differ  depending  on  the  material  constitutive  model. 


2  Kinematics  of  the  Single  Crystal 


The  basic  kinematics  of  the  representation  of  inelastic  behavior  under  finite  strains 
and  rotations  for  a  single  crystal  are  presented.  In  this  section,  two  “elastic” -“plastic” 
multiplicative  decompositions  of  the  deformation  gradient  will  be  discussed.  The  first 
representation  (referred  to  as  R-I)  obtains  the  relaxed  configuration  by  elastically  un¬ 
loading  to  a  stress  free  state  while  the  crystal  lattice  reorients  itself  with  its  initial 
reference  configuration.  In  R-I,  the  relaxed  configuration  is,  in  effect,  prescribed  by 
constitutive  constraints  on  the  plastic  spin.  The  second  representation  ( R-II)  obtains 
a  relaxed  configuration  by  elastically  unloading  to  a  stress  free  state  without  any  rota¬ 
tion.  This  requires  the  “elastic”  deformation  gradient  to  be  symmetric,  i.e.  a  stretch 
tensor.  The  spin  of  the  relaxed  configuration  of  R-II  will  be  shown  to  be  algebraically 
prescribed  due  to  the  symmetry  restriction  on  its  “elastic”  deformation  gradient.  These 
representations  are  discussed  in  detail  below  for  the  specific  case  of  the  single  crystal. 
Terms  which  are  specific  to  R-/will  be  subscripted  /;  terms  which  are  specific  to  the  R-II 
will  be  subscripted  II.  Terms  which  are  common  to  both,  such  as  the  total  deformation 
gradient,  F,  and  the  velocity  gradient,  L,  will  not  be  subscripted. 

Some  basic  kinematic  definitions  and  terminology  are  now  given.  As  discussed 
in  the  introduction,  a  material  point  in  its  original  undeformed  configuration  will  be 
represented  by  X.  The  same  material  point  in  the  deformed  configuration  will  be 
represented  by  x.  The  deformed  body  may  be  described  by  its  deformation  gradient, 
F,  which  is  given  by: 

F  =  Vx*.  (!) 

The  corresponding  rate  kinematics  begin  with  the  velocity  gradient,  L,  which  is  given 
by: 

L  =  Vxv  =  FF"1  =  D  +  W;  (2) 

where  D,  the  rate  of  deformation,  is  the  symmetric  part  of  L;  and  W,  the  spin,  is  the 
skew  part  of  L. 

2.1  Basic  Kinematics  of  Representation  I 

The  kinematics  which  are  specific  to  Representation  I  are  now  given,  beginning  with 
the  multiplicative  decomposition  of  the  deformation  gradient!  AS ARO  (1983)]: 
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F  =  FJFJ; 
F;  =  FJF/ ; 


(3) 

(4) 


where  F;  is  the  deformation  gradient  due  to  plastic  slip;  F*  is  the  residual  rotation 
of  the  lattice;  and  FJ  is  the  elastic  deformation  gradient.  In  other  terms,  the  gradient 
F}  accounts  for  the  stretching  and  rotation  of  the  lattice,  while  F/  accounts  for  plastic 
shearing  on  the  crystal  slip  planes  of  fixed  (reference)  orientation. 

The  rate  kinematics  may  also  be  decomposed  into  “elastic”  and  “plastic”  com¬ 
ponents.  This  is  done  by  substituting  the  decomposition  given  in  equation  (3)  into 
equation  (2),  from  which  we  obtain: 

l  =  fjf;-1  +  fjfjfj-'f;-1; 

Vj  =  f;f5f5_1fj"1  =  d?  +  wj;  u 

where  L*  is  the  velocity  gradient  of  the  relaxed  configuration  convected  to  the  deformed 
configuration,  F;  Dj  and  W*  are  the  rate  of  plastic  deformation  and  the  plastic  spin, 
respectively.  These  last  two  terms,  Dj  and  W*,  must  both  be  constitutively  prescribed 
in  this  formulation,  which  is  precisely  what  is  done  immediately  below  for  the  case  of 
the  single  crystal. 


Constitutive  Law  for  D/  and  Wj 

A  crystal  slip  system  may  be  defined  by  the  dyadic  s“  ®  m“,  where  the  vector  s“  lies 
along  the  a  slip  plane  and  the  vector  ma  is  perpendicular  to  this  slip  plane.  The  velocity 
gradient  in  the  relaxed  configuration  of  representation  I  is  given  by  Lj  =  F/F/-1.  This 
is  constitutively  obtained  by  the  summation  of  the  rate  of  plastic  straining,  -7“,  of  each 
slip  system  (s“  ®  ma): 

LJ  =  5^Tra*a®ma.  (6) 

a 

The  plastic  velocity  gradient  in  the  loaded  configuration  is  obtained  by  convecting  the 
slip  system  to  this  configuration: 

■**  ®  m*a  =  FJ(s“  ®  ma)FJ_1;  (7) 

resulting  in  the  following  expression  for  Lj  of  the  loaded  state: 

L?  =  £'yva®“*“-  (8) 
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The  rate  of  plastic  deformation  and  the  plastic  spin  in  the  loaded  and  unloaded  configu¬ 
rations  may  be  obtained  by  first  decomposing  the  plastic  slip  direction  tensors  sa  ®  m“ 
and  s*“  ®  m*®  into  their  symmetric  and  antisymmetric  parts,  which  we  define  by: 

S“  =  |(s®®m®  +  m“®8®); 

A*  =  -  m"®.°); 

S®  =  1(8*“®!!!*“  +  m*a®8*“);  1  ' 

A®  =  |(s*®  ®  m*®  -  m*®  ®  8*“); 

giving: 


D?  =  E„  T“S°; 
wj  =  Ea  TfaA4; 
D?  =  Ea  r s®; 

WJ  =  Ea  7“A®. 


(10) 


The  rates  Dj  and  W{  have  now  been  constitutively  prescribed  in  terms  of  the  plastic 
strain  rates  of  the  individual  slip  systems.  Later  in  this  paper,  a  specific  problem  con¬ 
sidering  a  “planar”  single  crystal  with  two  slip  systems  and  a  viscoplastic  constitutive 
law  for  Tf®  will  be  presented.  For  the  present,  we  can  see  that  this  first  representation 
accounts  for  the  material  and  lattice  deformation  of  the  crystal. 


2.2  Basic  Kinematics  of  Representation  II 


The  kinematics  specific  to  Representation  II  will  now  be  given.  The  multiplicative 
decomposition  of  RII  is  given  by: 

F  =  F^F',;  (11) 

F‘,I  =  F*„;  (12) 

where  F 7/  is  the  symmetric  elastic  deformation  gradient  (an  elastic  stretch  tensor); 
and  F7/  is  the  deformation  gradient  of  the  relaxed  configuration  obtained  by  elastically 
unloading  without  rotation  to  a  stress  free  state.  In  comparison  with  F;  which  describes 
the  plastic  shearing  of  the  material,  F7/  describes  the  plastic  shearing  as  well  as  any 
(residual  or  elastic)  lattice  rotations. 

The  velocity  gradient  is  decomposed  via  incorporation  of  equation  (11)  into  equation 
(2),  giving: 
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(13) 


L  —  IF*  IF*-1  _i_  TF«  T?P  TFP_11F*-1. 

”  *  II*  II  +  *  II*  II *  //  *11  t 

v„  =  f;,p57‘  =  dj,+  w;,; 

where  Lj,  is  the  velocity  gradient  of  the  relaxed  configuration  of  R-II ,  F,,;  D J,  is  the  rate 
of  deformation  of  the  relaxed  configuration  which  must  be  constitutively  prescribed; 
and  W J,  is  the  spin  of  the  relaxed  configuration.  The  spin,  Wj,,  is  algebraically 
prescribed  due  to  the  symmetry  restriction  on  F},.  This  is  shown  immediately  below: 

F}(  =  (D  +  W)FJ,  -  f;,(dj,  +  w;,), 
and, 

FJ,  =  (F;,)r  =  F}/(D  -  W)  -  (DJ,  -  W5,)F>„;  (14) 

(D  +  D5,)F'„  -  f;,(d  +  DJ,)  =  F5,(WJ,  -  w)  +  (w;,  -  w)f;,: 

therefore,  W*,  is  algebraically  obtained  to  be: 

WJ/-W-W//ID  +  DM;  (15) 

where  "Wu  is  a  fourth  order  tensor  with  components  of  order  elastic  strain  ONAT  (1987) 
mapping  symmetric  second  order  tensors  to  skew  tensors.  Therefore,  must  not  be 
constitutively  prescribed. 

Thus  far,  in  the  context  of  the  second  representation,  no  explicit  provision  has  been 
made  for  monitoring  the  orientation  of  the  crystal  lattice.  Onat  has  shown  theoretically 
that  R-II  can  accommodate  the  internal  state  and  orientation  of  a  material  element  by 
the  updating  of  an  appropriate  set  of  irreducible  even  rank  tensor  state  variables  ONAT 
(1982).  Here,  we  will  carry  out  this  proposal  for  the  problem  of  the  single  crystal. 

In  R-II,  the  deformation  gradient  of  the  relaxed  configuration  contains  both  the 
material  and  lattice  rotations  at  a  material  point.  Therefore,  for  the  case  of  the  single 
crystal,  F //  may  be  decomposed  into  the  product  of  lattice  rotation  RL  and  plastic 
slip,  FMp: 

FJ,  =  (16) 

Here,  F***’  describes  the  deformation  due  to  shearing  on  slip  planes  in  the  reference 
orientation  and  is  equivalent  to  FJ  of  R-I.  The  crystal  slip  system  in  the  relaxed  con¬ 
figuration  is  given  by: 

•*,  ®  m =  R^s"  ®  R1ma;  (17) 
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where  s“  and  ma  describe  the  initial  directions  of  each  slip  plane  and  its  normal.  The 
slip  system  may  be  further  decomposed  into  symmetric  and  skew  parts: 

S n  =  \  [s?/  ®  man  +  mau  ®  s?r] ;  (18) 

A®/  =  -[s“/  ®  m//  —  m//  ®  *//!• 

The  tensors  S7/  and  A JT  are  the  state  variables  needed  in  R-IIto  appropriately  monitor 
the  lattice  orientation.  It  is  emphasized  that  Sa7  and  Aa7  provide  an  irreducible  set 
of  even  ranked  tensors  which  monitor  the  lattice,  where  S“7  is  a  symmetric  traceless 
2nd  order  tensor  and  A f7  is  a  2nd  order  skew  tensor.  The  rate  of  change  of  these  state 
variables  is  found  to  be  given  by: 

Sf7  =  [Ri(Ri)T]  Sfj  +  S?7  [Ri(Ri)r]r;  (20) 

A//  =  [RL(RL)T]  Af/  +  A®/  (Rt(RL)r]r;  (21) 

where  RX'(RI')T  is  the  lattice  spin. 

The  velocity  gradient  of  the  relaxed  configuration  II,  Lpn,  may  be  expanded,  using 
equation  (16),  as: 

Lpu  =  F//F/71  =  Rt(Ri)T  +  RiF*<v(F*,’,’)_1(Ri)r,  (22) 

where  we  know  the  velocity  gradient  due  to  plastic  slipping  is  constitutively  prescribed 
by  the  rate  of  shearing  on  the  slip  systems.  Therefore,  using  (22)  and  (18, 19),  the  rate 
of  deformation  and  the  spin  of  relaxed  configuration  H,  can  be  expressed  as: 

d J,  =  Ei*sf/i  (23) 

a 

and 

Wpn  =  R*(R*)r  +  £  tT  A?,.  (24) 

a 

Since  W77  is  algebraically  defined  via  equation  (15),  we  obtain  an  expression  for  the 
lattice  spin  from  equation  (24): 

nL  £  Ri(Rt)T  =  Wj,  -  A  ft.  (25) 
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Equations  (20),  (21),  and  (25)  may  then  be  used  to  update  the  state.  The  state 
variables  Sn  and  An  may  be  updated  via  an  incremental  orthogonal  transformation, 
Q,  corresponding  to  ClLAt,  where  the  corotational  rates  of  Sq  and  An  are  given  by: 


s  j,°  =  s?7  -  +  s?7nL  =  o,  (26) 

a?, D  =  Af,  -  n‘Af,  +  A?,nL  =  0,  (27) 

and  the  “intermediate”  rotation,  Q,  must  therefore  satisfy: 

Q  =  nLQ,  (28) 

QW  =  I-  (29) 

The  state  variables  S“7  and  A 7/  are  then  updated  via: 

(S//)«+AI  —  (30) 

(a//)‘+a«  =  Q«+ai(A//)iQ^.a<.  (31) 


The  preceding  equations  are  easier  to  visualize  if  we  simplify  to  the  case  of  no  elastic 
stretching.  In  this  case,  the  expressions  in  terms  of  the  first  representation,  where  FJ 
is  now  simply  a  rotation  R*,  become: 

F  =  R*FJ; 

L  =  R*R*T  +  R'F7F7_1R*r; 

D  =  R*D?R*r;  [  ' 

W  —  R*R*r  -f  R*  W7R*t. 

The  expressions  in  terms  of  the  second  representation,  where  F7/  =  I,  become: 

L  —  L};; 

D  =  Dj/  =  E«r  S?7;  (33) 

W  =  WJ,; 

n4  =  W-E.7aA  Jr 

For  completeness,  the  kinematic  relationship  between  representations  I  and  II  is 
now  given.  We  first  pictorially  examine  a  body  as  it  goes  from  its  original  undeformed 
configuration  to  a  deformed  configuration  as  shown  in  Figure  1.  Here,  the  “relaxed” 
states  as  defined  by  F7/  and  F7  are  also  depicted.  Relaxed  configuration  I  is  a  stress  free 
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state  containing  no  lattice  deformation  or  changes  in  orientation.  Relaxed  configuration 
II  is  the  stress  free  state  which  contains  any  changes  in  lattice  orientation. 

Once  again,  we  examine  the  elastic-plastic  decomposition  of  the  deformation  gra¬ 
dient  as  implemented  in  R-I,  F  =  F}F7.  The  lattice  deformation  gradient,  FJ,  may 
be  broken  down  into  its  left  stretch  and  rotation  tensors  via  the  polar  decomposition 
theorem  to  give: 


F;  =  V*R*; 

(34) 

F  =  V*R*F7; 

(35) 

where  V*  describes  the  elastic  stretch  of  the  material.  The  decomposition  F  =  F7/F7/, 
where  F'/  is  a  symmetric  stretch  tensor,  yields  a  unique  F*7.  Therefore,  since  both  V* 
and  F \i  describe  the  elastic  stretch  incurred  by  the  solid,  we  must  conclude  that: 

V-  =  F‘„.  (36) 

In  other  words,  if  we  elastically  unload  to  a  stress-free  configuration  via  F,_1  =  FJJ1, 
then  FJ7  is  unique.  Therefore,  is  we  elastically  unload  to  a  stress-free  state  via  V*-1,  it 
must  be  the  same  stress-free  state  described  by  F7/.  The  rotation  matrix  R*r  retains 
an  unloaded  stress-free  configuration  and  simply  rotates  the  lattice  to  its  appropriate 
stress-free  configuration  which  would  be  its  initial  undeformed  configuration.  Therefore, 


the  following  relation  is  obtained: 

F?,  =  R*F7;  (37) 

or,  alternatively, 

FJ  =  R’rFj,;  (38) 

where  R*  is  the  rotation  tensor  which  defines  the  lattice  orientation.  Therefore,  it  is 
clear  that  R*  and  R£  are  identical  tensors.  The  relationship  between  the  rate  quantities 
of  R-I  and  R-II  is  also  established: 

Vjj  =  nt  +  RiL$Rir;  (39) 

Wj,  =  nL  +  RlW'RlT;  (40) 

Dff  =  RLD$RlT;  (41) 
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A  physical  interpretation  of  the  above  equation  is  that  the  R^L/R1,7  rotates  the  ve¬ 
locity  gradient  of  the  R-I  unloaded  configuration  into  the  R-II  unloaded  configuration. 
The  accounts  for  the  rate  of  this  rotation  and  is  the  spin  of  the  lattice.  In  other 
words,  R1  is  effectively  a  time-dependent  observer  transformation  relating  two  different 
elastically  unloaded  configurations. 


2.3  Example  Problem:  Planar  Single  Crystal 


The  problem  of  the  planar  single  crystal  with  two  shear  slip  systems  and  no  elastic 
stretching  is  now  examined.  This  problem  has  been  analyzed  by  PIERCE,  ASARO,  AND 
NEEDLEMAN  (1983).  For  completeness  and  for  comparative  purposes,  their  solution 
will  be  repeated  here.  A  solution  for  this  same  problem  using  the  second  representation 
will  also  be  presented. 

We  begin  with  some  definitions  of  basic  parameters  to  be  used  in  the  problem  as 
illustrated  in  Figure  2.  The  primary  and  conjugate  slip  systems  are  defined  by  the 
vector  pairs  (s^m1’)  and  (sc,  me),  respectively.  The  position  of  the  primary  system  is 
given  by  the  angle  0  clockwise  from  the  axis  of  imposed  tensile  deformation,  ej.  The 
position  of  the  conjugate  system  is  given  by  the  angle  0  counterclockwise  from  the  axis 
of  imposed  tensile  deformation,  es.  Initial  conditions  on  these  angles  were  taken  to  be: 
0(0)  =  40"  and  0(0)  =  20°.  The  angle  0  in  the  figure  describes  the  orientation  of  the 
lattice  with  initial  condition  of  /?(0)  =  0.  Therefore,  the  lattice  orientation  tensor  R* 
is  found  to  be: 


cos  (3  —sin0 
sin0  cos0 


(42) 


The  rate  of  plastic  shear  straining  on  the  primary  and  conjugate  slip  systems  is  given 
by  ip  and  7*,  respectively.  The  visco-plastic  constitutive  law  for  these  7“  is  given  by 
the  power-law  relation: 


where  d°  is  a  parameter  such  that  7“  =  d“  when  r“  =  ga;  ra  is  the  resolved  shear 
stress  on  that  slip  system;  ga  represents  the  internal  structure  of  the  material  and  may 
evolve  with  strain  hardening,  whereupon  the  ga  become  internal  state  variables  which 
monitor  the  hardness  of  each  slip  system;  and  m  indicates  the  rate  sensitivity  of  the 
material.  The  internal  structure  is  taken  to  evolve  according  to: 
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(44) 


9a  =  £  hafli*] 

0 

where  hafi  is  a  hardness  slope  matrix  and  therefore  accounts  for  latent  as  well  as  self 
strain  hardening.  Following  PIERCE,  ET  AL.,  for  this  example,  we  take: 

hpc  =  hgp  —  <jhee ;  (45I 

hpp  =  h]  hcc  —  h\ 


where, 


h  =  (*°^_+roY))  ;  (46) 

where,  for  aluminum  2.8  wt.%  copper  alloys,  ho  =  8.9r0  and  r,  =  1.8r0  and  1.0  <  q  <  1.4 
Pierce,  et  al.. 

2.3.1  Solution  Using  Representation  I 

The  relevant  material  law  information  has  now  been  given  and  we  now  move  onto  the 
solution  of  pulling  this  single  crystal  in  tension  along  the  zj  axis.  Since  the  crystal  is 
subjected  to  simple  tension,  points  which  lie  along  the  zs  axis  will  remain  along  this 
axis  giving: 

Pe,  =  Aej  (47) 

where  A  is  the  stretch  in  this  direction.  Using  representation  I  results  in: 

tanp  =  %/%  (48) 


a—  la  Fr»  Fi\tFin 

%'lw 

> 

where  Fj  is  obtained  from: 


FJ  =  ®  m’  4-  qrV  ®  me]  FJ. 

These  equations,  together  with  the  resolved  shear  stresses 

rp  =  fain  (2* -20); 

Te  =  fain  (2V>  +  20) ; 
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are  used  to  update  the  state  of  the  material.  The  orientation  of  the  lattice,  /?,  as 
a  function  of  the  amount  of  plastic  strain  on  the  primary  system  was  found  for  the 
following  three  cases  (all  with  d  =  10-1aee-1):  m  =  50,  q  =  1.4;  m  =  50,  q  —  1.0; 
m  =  10,  q  =  1.0.  The  results  are  shown  in  Figure  3  and  were  originally  calculated  by 
PIERCE,  ET  AL..  We  note  here  that  0  is  obtained  from  F?  by  equation  (48). 


2.3.2  Solution  Using  Representation  II 

Kinematic  representation  II  will  now  be  used  to  analyze  the  “planar”  single  crystal, 
where  the  lattice  orientation  is  included  as  a  state  variable  by  integrating  equations 
(20)  and  (21)  as  suggested  in  section  2.2. 

Since  we  are  considering  an  elastically  rigid  single  crystal,  FJj  =  I,  and  the  plastic 


deformation  gradient  is  equal  to  the  total  deformation  gradient: 

FJ,  =  F.  (52) 

Similarly,  the  plastic  spin  is  equal  to  the  total  spin: 

W'j  =  W  =  «£eti/[FF-1J  =  s/fceti/fFJjFjf1].  (53) 

This  spin  does  not  account  for  the  lattice  orientation.  The  lattice  spin  is  given  by: 

nL  =  (wf,  -  ya'„  -  .  (54) 

The  lattice  orientation  may  then  be  updated  with  equations  (30)  and  (31)  with  the 
initial  conditions: 

(S?,),^  =  "(s'*  +  »“],  (55) 

(A fjjtma  =  ^[s*  ®  m*  -  m"  ®  ■“],  (56) 


Since  this  is  a  planar  problem,  the  lattice  orientation  may  also  be  found  by  simply 
updating  the  angle  0,  where  0  =  fltn.  These  equations  together  with  the  constitutive 
laws  (40  -  43)  stated  earlier  are  used  to  update  the  state  of  the  material.  The  solution 
for  the  lattice  orientation,  0,  as  a  function  of  the  plastic  strain  on  the  primary  slip 
system,  ~f,  was  found  for  the  same  three  cases  as  the  earlier  solution.  The  results  are 
plotted  along  with  the  earlier  results  in  Figure  3,  and  with  which  they  are  identical. 
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It  is  important  to  note  that  the  solution  processes  using  the  two  different  represen¬ 
tations  are  very  similar.  However,  they  are  operationally  different  in  the  manner  in 
which  the  lattice  orientation  is  taken  into  account.  The  kinematics  of  the  elastic-plastic 
decomposition  of  the  deformation  gradient  in  representation  I  specifically  includes  the 
lattice  orientation.  The  second  kinematical  decomposition  is  more  general  and  does  not 
explicitly  contain  the  lattice  orientation,  which  must  be  implicitly  retained  and  “con- 
stitutively”  updated  as  a  tensor  state  variable.  Therefore,  this  representation  requires 
an  additional  tensor  to  properly  describe  the  material  behavior  of  a  single  crystal.  In 
R-I,  a  description  of  the  lattice  orientation  is  not  given,  ».e.  knowledge  of  the  plastic 
deformation  gradient  is  needed.  However,  in  R-II,  the  complete  state  of  the  material  at 
any  point  in  time  is  completely  described  by  the  instantaneous  hardnesses  of  the  slip 
systems  and  the  lattice  orientation.  In  the  next  section,  the  case  of  the  glassy  polymer 
is  examined  in  the  context  of  different  kinematical  decompositions. 

3  Kinematics  of  the  Glassy  Polymer 

Below,  we  will  give  a  brief  outline  of  the  constitutive  model  for  the  large  inelastic 
deformation  behavior  of  glassy  polymers  proposed  in  PARKS,  ET  AL.  and  BOYCE, 
ET  AL..  The  reference  configuration  of  the  glassy  polymer  is  an  isotropic  state  of 
the  material  which  consists  of  randomly  oriented  chains  (Figure  4).  We  assume  that, 
during  plastic  flow,  these  molecules  deform  affinely.  The  relaxed  configuration,  Fp,  is 
obtained  by  elastically  unloading  to  a  stress  free  state,  and  physically  indicates  the 
degree  of  permanent  molecular  orientation  existing  in  the  material.  In  general,  we  can 
express  the  total  deformation  gradient  as  the  product  of  the  elastic  stretch,  a  rotation, 
and  the  plastic  stretch: 

F  =  V'RU'.  (57) 

The  rotation  tensor  R  may  be  decomposed  into  the  product  of  elastic  and  plastic  ten¬ 
sors,  R  =  R«R',  where  F*  =  V'R'  =  R‘U*,  and  F *  =  R'U'  =  V'R',  via  the  polar 
decomposition  theorem.  However,  the  “elasticity”  and/or  “plasticity”  of  the  rotation 
tensor  is  indeterminate.  If  R  =  RF  and  R*  =  I  is  chosen,  all  such  rotation  effects 
are  lumped  into  the  affine  plastic  deformation  response  of  the  material.  This  results 
in  a  symmetric  and,  therefore,  unique  elastic  deformation  gradient,  F'T  =  F*.  This 
formulation  is  kinematic  representation  H  discussed  in  section  2  above.  Alternatively, 


if  R  =  R*  and  Rp  =  I  is  chosen,  F*  is  again  unique,  but  not  symmetric.  A  third 
alternative  for  rendering  this  factorization  unique  is  to  obtain  the  elastic-plastic  de¬ 
composition  of  the  rotation  tensor  by  imposing  a  constraint  on  the  spin  of  the  relaxed 
state.  Indeed,  this  was  precisely  the  procedure  followed  in  representation  I  above  for  the 
single  crystal,  where  the  plastic  spin  was  constitutively  prescribed,  thereby  generating 
a  particular  relaxed  configuration. 

In  the  remainder  of  the  discussion  on  the  glassy  polymer,  the  modelling  of  the 
material  behavior  will  be  considered  using  two  of  the  above  kinematical  representations: 
1)  a  relaxed  state  which  contains  all  rotations  is  chosen,  t.e.  Rp  =  R;  2)  a  relaxed  state 
is  chosen  via  a  particular  restriction  on  the  plastic  spin.  The  first  choice,  R-II,  was  used 
in  PARKS,  ET  AL.  and  BOYCE,  ET  AL.  In  this  case,  the  plastic  spin  is  algebraically 
defined  as  discussed  in  section  2.  The  alternative  of  defining  a  relaxed  configuration  by 
imposing  a  constraint  on  the  plastic  spin  will  be  made  clearer  after  further  discussion 
of  the  material  model  using  the  kinematic  representation  where  Rp  =  R. 

3.1  Choice  Is  Rp  =  R 

For  the  case  where  we  select  Rp  =  R,  where  R  is  defined  in  equation  (57),  and  F*  =  F*T, 
it  was  shown  earlier  that  the  plastic  spin  is  algebraically  defined  via  equation  (15).  A 
rate  of  plastic  shape  change  must  be  constitutively  prescribed  for  the  material.  We 
note  that  Dp7  describes  the  rate  of  shape  change  of  the  relaxed  configuration  given  by 
the  kinematics  F7/.  The  magnitude  of  D/j  is  given  by  the  plastic  shear  strain  rate, 
7 p,  and  the  tensor  direction  of  D/j  is  specified  by  N  =  ^ T*',  the  normalized  driving 
stress  state,  T*,  at  a  material  point.  The  driving  stress  state  is  given  by: 

T*  =  T  -  jFJjBFJ^,  (58) 

where:  T  is  the  Cauchy  stress  tensor  given  by  T  =  j£*[/nV*j  [ANAND  (1979)],  £'  is 
the  fourth  order  isotropic  elasticity  tensor,  and  V*  is  the  left  elastic  stretch  tensor;  J 
is  the  volume  change  given  by  deiF'u\  and  B  is  the  back  stress  tensor  resulting  from 
strain  hardening.  The  back  stress  is  uniquely  related  to  the  inelastic  distortion  in  the 
polymer,  where  B  is  coaxial  with  the  left  plastic  stretch  tensor: 


Tff  _  VP  PP 
*  II  ~  V  Ill'll 

polar  decomposition , 

(59) 

VJ,  =  Off  i'% 

eigenvalue  decomposition, 

(60) 
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Bi  =  Bi(  A?), 


(61) 


where  Bi  are  the  principal  values  of  B  defined  in  BOYCE,  ET  AL.,  and  Af  are  the 
principal  values  of  Ap.  The  coaxiality  of  B  and  Vf,  is  a  key  result  of  the  assumption 
that  the  molecules  deform  affinely  during  plastic  flow.  The  effective  equivalent  shear 
stress,  r,  is  given  by: 

r  =  (62) 

We  point  out  that  this  is  essentially  a  kinematic-hardening  law.  The  plastic  shear  strain 
rate  is  constitutively  prescribed  and  may  be  functionally  given  by: 

Y  =  Y  (r,s,p,9) ,  (63) 

where  s  is  the  scalar  athermal  shear  resistance  of  the  material,  which  evolves  with 
plastic  straining;  p  =  —  £trT  is  the  pressure;  and  9  is  the  absolute  temperature.  A 
detailed  expression  for  is  given  in  BOYCE,  ET  AL.. 

This  has  been  a  very  brief  description  of  the  kinematic  hardening  model  for  the  large 
inelastic  deformation  of  glassy  polymers  proposed  in  PARKS,  ET  AL.  and  BOYCE,  ET 
AL.  The  assumption  of  affine  response  of  the  molecular  chains,  which  lead  to  the  back 
stress  coaxiality  with  the  plastic  stretch,  was  instrumental  in  selecting  the  kinematic 
representation  of  deformation.  When  the  driving  stress  state  is  obtained,  the  back 
stress,  which  is  computed  in  the  relaxed  configuration,  is  transformed  to  the  loaded 
configuration.  Therefore,  alternative  relaxed  configurations  which  contain  different 
degrees  of  rotations,  i.e.  Rf  /  R,  should  indeed  result  in  the  same  material  response, 
providing  the  alternative  measures  of  back  stress  thus  obtained  are  also  consistently 
transformed  to  the  loaded  configuration.  As  stated  earlier,  another  unique  choice  of 
kinematic  decomposition  is  obtained  by  taking  R*  =  R  and  Rp  =  I  and,  therefore, 
Fp  =  Up.  The  resulting  principal  values  of  the  back  stress  tensor  would  be  identical, 
but  the  tensor  would  be  coaxial  with  Up  and  would  have  to  be  transformed  via  R* 
when  computing  the  driving  stress  state  T*.  This  decomposition  would  also  sufficiently 
complicate  the  rate  kinematics  so  as  to  make  it  an  unreasonable  choice  for  this  material. 
A  more  reasonable  alternate  choice  is  to  impose  a  constraint  on  Wp,  which  is  discussed 
next. 
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3.2  Choice  2:  Prescribe  Wp 


We  may  also  define  the  kinematics  of  the  glassy  polymer  in  a  manner  similar  to  rep¬ 
resentation  I  of  the  single  crystal  where  the  relaxed  configuration,  F/,  was  chosen  via 
a  constraint  on  both  the  rate  of  plastic  deformation  D;  and  the  plastic  spin  Wj.  In 
that  case,  the  constraint  was  constitutive.  For  the  case  of  the  glassy  polymer,  where 
the  molecules  deform  affinely,  there  is  no  obviously  analogous  or  intuitive  choice  for  a 
“constitutive”  constraint  on  the  spin  of  the  relaxed  configuration.  Therefore,  we  may 
obtain  our  relaxed  state  represented  by  Fp  by  choosing  a  convenient  Wp.  Indeed,  note 
the  “convenience”  offered  in  R-I  by  the  imposition  of  the  particular  constraint  on  the 
spin  of  the  relaxed  configuration  of  the  single  crystal.  The  simplest  possible  form  for 
the  spin  of  the  relaxed  state  is  chosen: 

w;„  =  0.  (64) 

In  other  words,  a  relaxed  configuration  which  is  not  spinning  is  chosen.  In  general,  as 
in  R-I,  this  will  result  in  elastic  and  plastic  deformation  gradients  which  both  contain 
rotations,  F  =  FJ/7Fpr/  =  V*RJJJRP//UP.  The  rate  of  deformation  of  the  relaxed 
configuration  is  given  by: 

DJ„  =  I'M,  (65) 

where, 

M  =  -^=-RJ7J[T  -  jFJ/7fiFJ77]R*/7;  (66) 

where  6  is  coaxial  with  the  left  plastic  stretch  of  F ///,  V///  =  R///UpRj//t;  T  = 
i£*[/nV*];  and,  therefore,  M  is  the  direction  of  the  driving  stress  state  in  the  relaxed 
configuration,  F7//.  We  note  that  the  principal  values  of  B  are  identical  to  those  of  B 
given  above  in  equation  (61). 

A  schematic  of  the  deforming  glassy  polymer  is  shown  in  Figure  5.  Here,  the  relaxed 
configuration  of  both  representations  II  and  m  are  shown.  We  note  that  the  two  relaxed 
configurations  simply  differ  by  a  rotation,  where  we  have: 


P  J,  =  KU»; 

(67) 

K,  = 

(68) 

which  results  in: 
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FP  p«  TfP 

n  —  nr 


(69) 


The  corresponding  spin  of  these  two  re  ixec  'onfigurations  are  given  by: 

I*///  —  Djii  =  ^ni^ni  5  (70) 

and 

L//  =  D//  +  Wj,  =  RJ//RJ//  +  ®7//^//®7//-  (71) 

Therefore,  W7/  =  R///R///.  This  connection  between  the  two  relaxed  configurations 
for  the  glassy  polymer  is  analogous  to  that  presented  in  section  2  for  the  single  crystal. 
A  detailed  proof  of  the  independence  of  the  solution  upon  the  kinematic  decomposition 
due  to  the  isotropic  character  of  both  the  Cauchy  stress  and  the  back  stress  is  given  in 
the  Appendix. 

Identical  results  are,  of  course,  obtained  using  the  kinematic  decomposition  of  rep¬ 
resentations  II  and  in  for  the  glassy  polymer.  In  particular,  the  analysis  of  the  simple 
shear  of  a  glassy  polymer,  polymethylmethacrylate,  was  conducted.  The  functions  for 
ip  and  B  and  all  necessary  material  properties  are  given  in  BOYCE,  ET  AL..  The 
resulting  shear  stress-displacement  curves  for  both  R-II  and  R-III  are  shown  in  Figure 
6.  The  solutions  are  identical. 

4  Conclusion 

This  paper  has  examined  the  representation  of  the  kinematics  of  inelastic  deformation 
problems  involving  finite  strains  and  rotations.  The  independence  of  problem  solution 
on  kinematic  representation  was  demonstrated  using  two  distinct  examples.  The  first 
example  considered  the  planar  single  crystal  of  PIERCE,  ET  AL.  subjected  to  tensile 
loading.  The  second  example  analyzed  the  glassy  polymer  as  constitutively  modelled  by 
BOYCE,  ET  AL.  subjected  to  simple  shear  loading.  Both  problems  were  solved  utilizing 
various  representations  of  the  kinematics.  The  solutions  were  found  to  be  independent 
of  the  choice  of  the  relaxed  configuration,  i.e.  the  representation,  as  shown  in  Figures 
3  and  6.  The  two  cases  discussed  in  this  paper  were  examined  with  physically-based 
constitutive  models,  where  the  tensorial  representation  of  the  resistances  to  deformation 
in  the  single  crystal  and  the  glassy  polymer  are,  in  a  sense,  clearly  understood.  The 
tensorial  state  representation  of  the  physics  of  deformation  is  not,  at  this  point,  as 
well  characterized  for  more  complex  materials  such  as  polycrystalline  metals.  This  has 
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lead  to  the  corresponding  debate  over  the  representation  of  the  kinematics  of  plastic 
deformation  in  the  presence  of  anisotropy  and  the  role  of  the  “plastic  spin” . 

It  can  be  concluded  that  the  most  convenient  representation  for  the  material  of  con¬ 
cern  may  be  chosen  when  solving  problems  involving  finite  strain  plasticity.  However, 
on  reflection,  one  also  wishes  to  acknowledge  that  it  is  of  great  interest  to  retain  in¬ 
formation  on  the  state  and  orientation  of  the  material.  General  representations  which 
monitor  this  information  with  internal  state  variables  are  highly  desirable.  General 
representations  also  permit  a  greater  freedom  for  improvement  or  expansion  of  existing 
material  models  by  providing  a  kinematic  framework  independent  of  material  model 
which  will  not  require  modification  to  accomodate  material  law  modifications. 
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Figure  2.  Schematic  of  the  planar  single  crystal  containing  two  slip 
systems.  The  primary  and  conjugate  systems  are  depicted  with  respect 
to  the  tensile  load  axis. 
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Figure  3.  Lattice  orientation,  0,  vs  primary  plastic  shear,  yp, 
using  both  representations  I  and  II.  Results  are  identical. 


Figure  4.  Schematic  of  amorphous  polymer. 


Figure  6.  Simple  Shear  of  PMMA  —  Shear  Stress  vs.  Normalized 
Displacement  using  kinematic  representations  II  and  III.  Results 
are  identical. 


Appendix  A 


In  section  3  of  this  paper,  a  symmetric  form  for  the  plastic  velocity  gradient  was 
postulated.  The  resulting  decomposition  was  referred  to  as  representation  III,  R-III.  It 
is  the  purpose  of  this  appendix  to  show  that,  in  general,  the  introduction  of  an  arbitrary 
skew-symmetric  component  of  the  plastic  velocity  gradient  tensor  does  not  affect  the 
Cauchy  stress  material  response  as  long  as  the  remaining  constitutive  equations  are 
form  invariant  as  in  the  case  of  the  glassy  polymer.  In  the  paper,  such  an  example  was 
analyzed,  where  the  “Wp”  which  was  an  algebraic  result  of  the  symmetry  imposed  upon 
F*  was  the  “arbitrary"  skew-symmetric  component.  The  corresponding  decomposition 
was  referred  to  as  representation  II,  R-II. 

For  the  proof,  we  consider  two  constitutive  initial  value  problems,  where  the  only 
difference  is  the  flow  rule: 

Case  a)  The  plastic  velocity  gradient  as  expressed  in  the  relaxed  configuration  is  taken 
to  be  symmetric  (RIII  of  the  paper). 

L^sD^T.S,*); 

**///  = 

3  =  g(T,B,3);  (A.l) 

T  =  r  [InUJ„]; 

fi  = 

with  the  initial  conditions  of  F7//(0)  =  I,  (and,  therefore,  B(0)  =  0)  and  3(0)  =  s0- 
The  overlaid  bar  indicates  that  these  are  expressions  for  the  description  of  state  in  the 
relaxed  configuration,  F7//,  t.e.  not  in  the  loaded  configuration.  Note  that  Dp,  g,  £*, 
and  f  are  isotropic  functions. 

Case  b)The  plastic  velocity  gradient  is  taken  to  be  arbitrary.  (One  such  example  was 
R-II  given  in  the  body  of  this  paper.) 
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L»  =  D*(T,B,3)  +  Wj; 

Fj  =  [D*(T,B,3)  +  WjjFj; 
s  =  g(T,B,3); 

T  =  £*[ZnUJ]; 

B=f(Vf); 


(A.2) 


with  the  initial  conditions  of  F*(0)  =  I,  (and,  therefore,  B(0)  =  0)  and  3  =  a0-  The 
overlaid  tilde  indicates  that  these  are  expressions  for  the  description  of  state  in  the 
relaxed  configuration,  Fj,  ».e.  not  in  the  loaded  configuration. 

In  order  to  prove  our  assertion  that  the  system  of  equations  in  cases  (a)  and  (b) 
will  lead  to  identical  solutions,  we  consider  a  “trial”  solution  for  the  second  system  in 
the  form: 

fLl;QF"';  (A.3) 

where  Q  is  a  rotation  tensor  to  be  determined.  From  these  relations,  we  obtain: 

r;  =  R-,„qr; 

u;  =  QV,„Qt;  (a.4) 

vf  =  qv;„Qr. 

Due  to  the  isotropic  character  of  f  and  £*,  the  above  relations  yield: 

/n(UJ)  =  Q/n(U‘m)Qr; 

B  =  QBQr;  (A.5) 

T  =  QTQr. 

Next,  differentiating  Fj  as  given  in  equation  (A3),  and  then  substituting  in  the  flow 
rules  of  cases  (a)  and  (b),  we  obtain: 


(A.4) 


(A.5) 


F»FJ_1  =  QQr  +  QD'ff  ,B,  J)Qr  =  D»(T,B,S)  +  WJ. 


(A.6) 


Due  to  (A5)  and  the  symmetry  and  isotropic  character  of  D*,  equation  (A6)  is  satisfied 
provided: 


QQr  =  Wf. 

From  equations  (A4)  and  (A5),  we  also  verify  that: 
(<kiu;)-'R;TR;r  =  (*tu;„)-'R;„TR;,/. 


(A.7) 


(A.8) 


Therefore,  the  Cauchy  stress  tensor  T  obtained  by  integration  of  one  or  the  other 
system  of  constitutive  equations  is  identical. 
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